Galactic Halos are Einstein Clusters of WIMPs 
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It is shown that an Einstein cluster of WIMPs, WIMPs on stable circular geodesic orbits gener- 
ating the static spherically symmetric gravitational field of a galactic halo, can exactly reproduce 
the rotation curve of any galaxy simply by adjusting the local angular momentum distribution and 
consequent number distribution of the WIMPs. No new physics is involved (assuming the forth- 
coming discovery of WIMPs) . Further, stability of the orbits can require an inner truncation of the 
halo and an explicit example of this is given. There is no exact Newtonian counterpart to the model 
presented since pure Newtonian gravity fails to account for the contribution made by the angular 
momentum distribution of the particles in creating the gravitational field. In effect, a galactic halo 
is supported by the hoop stresses created by the orbiting WIMPs. 



Many years ago Einstein Q] showed that one can in- 
troduce rotation without global angular momentum. He 
used a cluster of randomly distributed gravitating dust 
(non-colliding) particles each on a circular geodesic orbit 
about the center of symmetry of the gravitational field 
produced by the entire distribution of particles. In this 
letter I show that this model, in a modern context, can 
solve all dynamical issues associated with galactic rota- 
tion curves, one of the foremost problems in all of physical 
science today. It is proposed that the dark matter halos 
are the gravitational fields produced by massive parti- 
cles (WIMPs) in orbit about galactic centers |2j. The 
gravitational contribution of the visible components is, 
in excellent approximation, simply ignored here. 

Einstein's heuristic arguments can be made more pre- 
cise 3] and here we simply sketch the construction of the 
energy-momentum tensor. Consider a static spherically 
symmetric continuous distribution of particles (WIMPs) 
of (rest) mass m/0 and 4 - momenta p a which are on 
stable circular geodesic orbits about the center of sym- 
metry. (The existence and stability of these orbits is dis- 
cussed below.) In terms of "curvature" coordinates the 
template for our calculations takes the form Q| 



ds 2 



dr 2 



2M(r) 



2 dn 2 ~ e 2 *^dt 2 , 



(1) 



where dQ 2 is the metric of a unit sphere (d8 2 +sin 2 
Now let n be the local number density of particles and 
at some spatial point V (labelled, say, by {r v ,9 v ,<j) v )) 
construct 



<p a p P > 



(2) 



where the average < > is taken at V and over all parti- 
cles passing through V. The total angular momentum L 
(defined as usual by 1? = p 2 , +p 2 / sin 2 9 ^ 0) is of course 
constant along each trajectory but we also assume that 
it is constant along every trajectory through V thus ar- 
riving at the main results 



< Pt, > - 



< P<j> >- 



P* >= o 



(3) 



and 



< vl >- 



<P%> 
sin 2 6„ 



L 2 
2 ' 



(4) 



Moving to a continuous distribution and maintaining 
spherical symmetry we have L = L(r) and n = n(r) 
so that finally 



T\ = -mn(r){l + C 2 (r)) 



and 



where 



Tt = Tl 



C(r) 



mn(r) r2 



L(r) 



(5) 



(6) 



(7) 



All other components of T@ vanish and V^T^ vanishes 
identically. 

Let us now make use of the standard Einstein equations 
Ga = ^ 7r ^a m the coordinates defined by (JTJ and using 
the energy-momentum tensor derived above. First, since 
TJT = 0, we can solve for M{r) algebraically from G r r = 
to obtain 



M(r) 



7^<J> 



1 + 2r$' 



where = d/dr. Next, writing ~TljT e e 
arrive at 



r$ = 



C 2 {r) 



(8) 

~G\/Gl we 



(9) 



l + £ 2 (r) 

so that ||SJ can be written in the equivalent form 

M (r) = 

V ' l + 3£ 2 (r) 

Finally, from JS| and © since T t * + 2Tg = —ran, we have 
(1 - r$')(r$" + 2$' + 2r($') 2 ) 



(10) 



Airmn 



r(l + 2r$') 2 



(11) 
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which, with the aide of ©, can also be given in the form 



4nmn — 



C(2rC + £(l + 3£ 2 )) 
r 2 (l + £ 2 )(l + 3£ 2 ) 2 ' 



(12) 



Of course, van has to be > 0. In this model there is no 
way to determine m and n separately (an advantage dis- 
cussed below). We refer to the product nm as the "num- 
ber distribution" . As l|12|) makes clear, in this model the 
number distribution mn can be considered a consequence 
of the angular momentum distribution C. 

If we were given C(r) then the gravitational field of a 
halo would be complete up to quadrature in 10. As dis- 
cussed in detail previously it is $(r) (and not C{r)) 
which follows directly from observations of galactic ro- 
tation curves. However, important constraints can be 
placed on C{r). In the notation of 0, the existence of 
these circular orbits requires < r$ < 1 which with © 
translates simply to < C(r) 2 < oo but shows that the 
first term in (|llfl is positive. The stability of the circular 
geodesic orbits requires 3$ +r$ > 2r($ ) 2 which gives 
rCC > ~ C 2 and shows that the second term in i|llf) is 
positive only for r<!> > 1/4 (C 2 > 1/3). A mapping onto 
the observer's plane requires <& > r<i> + 2r($ ) 2 so that 
rCC < C 2 and so the distribution of angular momentum 
should satisfy 



L 2 < rCC < C 2 



(13) 



which, when applied to l|12|l directly repeats the impor- 
tant result: The stability of the circular geodesies orbits 
requires 



C 2 >1 



(14) 



Let r signify the border of stable orbits defined by 
£ 2 (ro) = 1/3. If C 2 (r) is monotone increasing then the 
region r < ro must be excluded and the halo distribu- 
tion is necessarily a thick shell. An explicit example of 
this is shown below. (If one was to measure the "rel- 
ativistic" nature of the configurations by way of (say) 
2M(r)/r then they are always highly "relativistic" , start- 
ing with r ~ 6M(r) and becoming more "relativistic" 
(r ~ 3M(r)) as C 2 — > oo. This notion of "relativistic" 
is, however, based on the Schwarzschild vacuum (M = 
constant) and it is inappropriate here. This is discussed 
further below.) 

As discussed previously p| , observations of galactic ro- 
tation curves are reported by way of the "optical conven- 
tion" v = (A — A e )/A e so that 



$ b 



'7,2 



r(l - r$') 



(v(b) - v(b = 0)) 2 



(15) 



where the mapping from the observer's plane to r is given 
by way of the classical impact parameter 



r 

2<t> ' 



(16) 



Observations therefore give <& up to quadrature and as 
a consequence the full geometry, the angular momentum 
distribution C and consequently the number distribution 
mn. In any particular case one can take the view that the 
angular momentum distribution of the orbiting WIMPs 
has been adjusted to exactly match the observed rota- 
tion curve. Equivalently, the observed rotation curve is 
a consequence of the distribution C. 

With the understanding that the procedure used here 
can be applied to any rotation curve, we consider, as 
previously , the specific example of a universal curve 
given by 



v{bf 



ab 2 



b 2 +f3' 



(17) 



where a and (3 are positive constants characteristic of a 
particular galaxy. With (|17fl it follows form l|15l) that 



„2$ 



ar 2 W(y) - (3 
where W is the Lambert W function Q, 

8 



y 



ce > 



(18) 



(19) 



and c is a constant > 0. Observe that with (|18f) and H19|) 

1 



l + W 



(20) 



and so from (|5J the required angular momentum distri- 
bution is given simply by 



L 2 = 



1 



(21) 



The gravitational mass follows from (JHJ and is given by 



M = 



(22) 



3 + W 

Finally, from (|llfl or (1121) the number distribution follows 



W(ar 2 (3 + 6W + W 2 ) +2(3W) , , 

4nnm= ^ ^ ^ ^ — '-. (23 

(3 + W) 2 (l + W) Va 

In this model, mapping onto the observer's plane requires 
c > a(3. With (|21|) , limits on W follow from the stability 
condition (|14|) and the existence condition for the orbits 
so that 



< W < 3. 



(24) 



Since W(y) is a monotone decreasing function of r, de- 
creasing from oo at r = to as r — > oo, condition 124|) 
tells us that for stability the distribution is necessarily 
a thick "shell" with the central divergence in W(y) re- 
moved as mentioned above. The junction conditions as- 
sociated with this inner truncation are discussed in detail 
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FIG. 1: Some representative properties of the model which 
follows from 11711 . For this figure a = f3 = 1 and c — 2. These 
values are for demonstration purposes only. The thick solid 
curve gives W(y), and the thin solid curve gives M(r) which 
is essentially linear. The dashed increasing curve gives C? 
and the lower dashdot decreasing curve gives 87rmn(r) which 
lies below the effective density, the dashdot decreasing curve 
given by 8nmn(r)(l + C 2 ) (in these two cases 87r is included 
as an amplification factor for the image). The distribution is 
truncated when the orbits become unstable at VV = 3, in this 
case for r ~ 0.633. The curves labelled N and B are discussed 
below. 



below. Some representative properties of this particular 
model are shown in Figure ^ 

Because of the stability condition I] 14(1. a discussion 
of junction conditions is appropriate. Geometrically, the 
junction is examined by way of the Darmois-Israel con- 
ditions - the continuity of the first and second funda- 
mental forms at a boundary surface £ 8]. For example, 
the smooth junction of metrics of the form QJ, at any 
boundary surface defined by constant r, only requires 
the continuity of M and $ , assuming the continuity of 
r, and 4>. As a result, Gy, but not Gg nor G\, is nec- 
essarily continuous at rs- As a consequence, since any 
Einstein cluster has G r r = 0, the matching conditions be- 
come trivial. One could, for example, insert an interior or 
exterior Schwarzschild field at any simply by choosing 
the Schwarzschild mass M to be 



M 



1 + 2r$' 



(25) 



This "Swiss-cheese" type of matching is, however, clearly 
inappropriate for the objects considered here |9j]. Here 
the distribution of "dust" particles is supported by the 
hoop stresses given by (JBJ and junctions onto other forms 
of "dust" without these stresses are more appropriate 
in the present context. Consider, as a simple example, 
junction onto a uniform dust section (in particular, a 



piece of a flat Robertson- Walker spacetime with density 
p(t) but not necessarily the background cosmology). 
Let the junction surface £ be defined by Mj r = e where e 
is a (dimensionless) constant. The continuity of the areal 
radii and effective gravitational masses then gives 



r s = a(t)r s = C\j —a(t) 



\3/2 



(26) 



where c is the velocity of light and C is the constant 
defined by 47rp(<)a(i) 3 /3. Unlike the "Swiss-cheese" type 
of model, here both rs and r^ evolve with t. 

It is interesting to compare the model presented here 
with what would result from pure Newtonian gravity. To 
do this consider the potential $ and density p defined by 

[HI 



V 2 $ = Airp : 



(27) 



where R\ is the time component of the Ricci tensor of 
(JTJ. A direct computation shows that 



1 + r$ 
1 - r<5' 



(28) 



where mn is given by (|11J) . or equivalently from 

p = mn(l + 2C 2 ). (29) 



In pure Newtonian gravity we have 
p = mn 



(30) 



and so pure Newtonian gravity fails to account for the 
contribution made by the angular momentum distribu- 
tion of the orbiting particles. This contribution is in fact 
the sole support for the Einstein cluster. 

We end with a few words about gravitational lensing. 
Of course any model of a galactic halo should not only 
reproduce the observed rotation curves, but also agree 
with any observed lensing. At present the possibility 
of such simultaneous measurements seems remote |l2l. 
Here we only demonstrate that an Einstein cluster can be 
expected to be a significant source of gravitational lens- 
ing. The curves B and A*" in Figure ^show the potential 
impact parameters B [l3T | for non-radial null geodesies in 
the Newtonian case (N, no bending) and for the Einstein 
cluster described above (B, for the values of a, (3 and c 
stated). The flattening of B is clear from (I22|l . 

In summary, we have presented a model, based on an 
argument originally due to Einstein 1], that exactly re- 
produces the rotation curve of any galaxy. The model 
requires no new physics but does require the existence 
of WIMPs on stable circular geodesic orbits. It is these 
orbiting particles that create the gravitational field of 
the dark matter halo. The halo is held together by the 
hoop stress generated by the angular momentum of the 
particles. The actual value of the rest mass for these 
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particles plays no role in the argument which, given that 
this is not presently know, can be considered an advan- 
tage. There is no exact Newtonian counterpart to the 
model since pure Newtonian gravity fails to account for 
the contribution made by the angular momentum distri- 
bution of the particles in creating the gravitational field. 
This contribution is in effect the entire content of the 
model. The model predicts, on the basis of the stability 
of orbits, that the halo is shell-like and so the WIMPs do 
not occupy the inner regions, at least on stable circular 
orbits. The gravitational lensing properties of the model 
need to be compared with observations when these be- 



come available. These properties involve an analysis of 
the junction of the halos onto both interior and exterior 
dust fields which are not supported by hoop stresses. 
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